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1. INTRODUCTION 

Remotely operated vehicles (ROVs), like all underwater vehicles [1], have grown in importance 
over the years. Their deployment is becoming more and more omnipresent and often unavoidable [2], 
because the underwater environment is a hostile and formidable environment and makes this type of mission 
difficult and expensive. Marine robots provide the more than perfect solution for missions [3], involving civil 
and military fields such as maritime cartography, oceanography, underwater research or seabed surveillance. 

Several models and technologies of these machines were appeared, at their beginning, they were 
remote controlled with the humane presence in the chain of decision, after they created autonomous variants 
AUVs (autonomous underwater vehicles) which filled the limitations of ROVs. All these submersibles have 
been able to assert their importance and consolidate their position as an indispensable tool to avoid 
endangering humans and to reduce operating costs by automating some repetitive or reduced-decision tasks 
in relation to the fields aforementioned [4]. 

The propellant solution for torpedo-shaped, used in our work, underwater robots is commonly a rear 
thruster as shown in Figure 1, while the guiding solution is almost generally a combination of rudders and 
diving bars to ensure the evolution of the submarine in the aquatic environment in the desired directions. An 
architectural solution, as shown Figure 2, consists in using four orientable propellers with a pair at the front 
and another pair at the rear, but this solution supposes the installation of additional mechanisms to handle the 
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thrusters and in addition, it also has a Maneuvering thruster in bow. But is it possible to improve or optimize 
this architectural solution and offer a more advantageous alternative for shifting? 

The architectural solution proposed in this study consists of a torpedo-shaped marine robot, but with 
four propeller thrusters at the rear, which will simultaneously provide propulsion and direction guidance. The 
principle of this solution is partly inspired by aerial drones with four vertical rotors and at the same time also 
used as an emergency man oeuvre for surface navigation in the event of rudder failure, but only for ships 
with more than one propeller. So, to keep a given direction, it is necessary to maintain the same thruster 
speed and for changing direction it will only be necessary to unbalance the distribution of the thrust of the 
thrusters in favor of one direction or the other. 


Figure 1. Submarine robot with rear propeller, rudders, Figure 2. Alistar robot with four rear propellers 
and diving bars (ECA group) (ECA group) 


2. RESEARCH METHOD 
2.1. Presentation 

The general concept of this vehicle was presented in the article [6] which precedes this study, it 
describes in addition to the shape of this vehicle, the advantages and disadvantages of this architectural 
solution. In the same way, we were able to demonstrate in our last publication [7] that in a real aquatic 
environment, this ROV as Figure 3 can move and execute the orders received. Thus, as demonstrated, the 
most common engine thrust combinations have been tested and the resulting curves have been presented. 


Figure 3. Proposed model in perspective 


Two approaches are possible to have a dynamic model [8]: one starting from the conservation of 
equations of the kinetic and potential energies of Lagrange-Euler [9], [10] and the other developing the 
fundamental equation of mechanics named method of Newton-Euler [11]. It's this last that was used in this 
researched works. The modeling of the ROV includes the kinematic part that will focus on the movement and 
geometric relationships of the submarine. While the dynamic part will deal with the forces and torques acting 
on this machine. 


2.2. Kinematic model 
The global position vector [12], [13]: 


n = [x,y,z,9,9,]" (1) 


n = non)" n, = x,y,z)" n, = (p,8, p)" (2) 
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In the same way we can combine the speeds: 

v = wv)" v = (u v,w)” v = (p, qr)" (3) 
In coordinate system Ro: 

I, =Ja (mM f = Je, (1,)V2 (4) 
with: 


cpc0 —swWcpt+cpsOsp spsp +cwpsOcp 
(5) 


Ja (m) = (s0 c0 cicyt+sWsOsp — cy so + sy sh co 
— sé c0 sp c0 cp 


1 singtan@ cosgtané 


Jalm) = fo coso — sing (6) 
0 sing/cos0 cosg/cosé 


The Roll, Pitch and Yaw (9, andy) designates the Euler's angles (Figure 4), which in robotics 
correspond to the system commonly called R.T.L: 


Figure 4. Representation of Euler angles [14] 


2.3. Dynamic model 

In general, the forces generated by the rotation of the earth on the vehicle can be neglected 
compared to the hydrodynamic forces, if we consider that the Ro coordinate system is a Galilean coordinate 
system fixed to the surface of the Earth. The dynamic simulation will be inspired by the work of Fossen [9]. 
In this reference, the fundamental principle of the dynamics applied to the mobile gives [15]: 


Mv=f.+fotfrtrtt+w (7) 


Where the vectors: 

— v=[u,v,w,p,q,7]' : vehicle accelerations 

— fe: forces and torques of drive inertia and Coriolis 

— fg: forces and torques induced by weight and Archimedes' thrust 
— fn: hydrodynamic forces and torques 

— t: forces and torques produced by actuators 

— w: external disturbances (waves ...) 


2.3.1. Inertia and coriolis 

Any moving body moving in a fluid whose densities are close, causes a displacement of a certain 
quantity of this fluid during its movement. Thus, from the physical point of view, Newton's second law is no 
longer applicable in his classical form, because air was admitted as a vacuum and the mass of the displaced 
air is negligible compared to the mass of the solid in motion, from where it is necessary to take into account 
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this displaced mass of fluid [16]. Consequently, the balance of the forces due to the inertia and to this added 
mass of water must be written in the form: 


È fi = (M, + Ma). Yr (8) 


As a result, relation (8) which concerns the mass matrix, is seen to be developed. The inertia matrix 
of the fluid is positive and symmetrical. The underwater vehicle has two planes of symmetry, one with 
respect to the plane (xz) and another with respect to the plane (xy) and since it has a torpedo shape, so we can 
simplify the symmetries and the matrix can take the following form [17]: 


X, 0 0 0 0 0 
oY% 0 0 0 0 
0 0z% 00 0 
Ma=-|o 0o 0 K& 0 0 (9) 
0 0 0 0 M 0 
0 00 0 0 wẹ 
For our ROV: 
Xa=0 1m 2.86 
‘nee, 2 
se ie 88.84 
pe 2 
lap = praeh 88.84 
pral? ma’ 
45—22 TIEL 10.66 
273 3 
2 peek a 10.66 


E 12 15L 


Likewise, for the matrix of inertia and as the center of the reference, Rv coincides with the center of 
gravity and taking into account the symmetries we can simplify is being as [17]: 


m 0 0 0 0 0 
0 m0 0 0 0 
0 0m 0 0 0 

M =-j|o 0 0 ly 0 0 (10) 
G:F 200. 20 E 
0 0 0 0 0 Iz 


2.3.2. Weight and Archimedes' thrust 
The forces that the submarine undergoes in water are Archimedes' push and that of its own weight 
and are written: 


W=m.g (11) 
B = p.V.g (12) 


with: 
— m: vehicle mass 
— g: earth acceleration 
— p: water density 
— V: displaced water volume. 
Thus, the vector of hydrostatic forces can be written as [16]: 
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—(W — B) sind 
(W — B) cos@ sing 
(W — B) cos@ coso 
fs = (Yg W — YpB) cos@ cosp — (zW — ZB) cos® sing 
—(x,W — x,B) cos@ cosp — (zgW — z,B) sind 
(xgW — x,B) cos@ cosp — (YW — ypB) sind 


(13) 


with: 
— B: Archimedes' thrust (Buoyancy) in N 
— W: Weight in N 
— (%g,Xg,Xq): the position of the center of gravity in Rv 
— (Xp, Xp, Xp): the position of the center of gravity in Rv 
If we consider that the center of thrust and the center of the reference Rv, are placed on the same 
axis (Oz), then, the center of gravity must be below the center of thrust for the submarine to keep its position 
of initial equilibrium under the effect of the righting moment (Figure 5). 


Figure 5. “x” configuration in rear view 


From where x, = 0,y, = 0 and Zp = —h,and that the mass of the submarine is distributed 
symmetrically with respect to the three planes (xz), (xy) and (yz), we know that in Ry. xg = 0, yg = 
Oet Zg = 0, the same at equilibrium W=B: 


0 
0 
0 
W h cos® sing 
W hsin@ 
0 


fa = (14) 


2.3.3. Hydrodynamic forces and torques 

There are other types of hydrodynamic dampers and they affect primarily surface ships rather than 
underwater vehicles. Thus, it becomes difficult to separate all dampers closely, it is then necessary to express 
in a global way all these hydrodynamic damping forces in a single term which depends on the speed [18]: 


fn = —D(v).v = —(D, + Dy (v)).v (15) 
— Dy: linear damping. 


— Dyw: non-linear damping. 
with: 
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D 0 0 0 0 0 
0D 0 0 0 0 
0 0 Dy 0 0 0 
Pe) 0% o 0 pO 0 (16) 
0 0 0 0 D 0 
0 0 0 0 0 D 
Dyulul 0 0 0 0 0 
O Dwbwl 0 0 0 0 
0  Dwlwi 0 0 0 
Du=| 9 0 i Daplpl 0 0 a7) 
0 0 0 0 Dnglal 0 
0 0 0 0 0 Darlt| 


Where the coefficients in « | | » are the absolute values of the linear and absolute speeds. 
The submarine will be assumed like a cylindrical shape of length L, with two hemispheres at the end 
of radius a. As the submarine does not have a drive, rudder or rudders therefore: 


D, == = pr(2aL + na?) Drult| = 5 pra? (18) 
D, === pra? Dnlv| = pal (19) 
Dy, =0 Daw|w| = pal (20) 
Dy = 0 Dnp|P| = $ pal (21) 
D, = 0 Dnqlal = = pal’ (22) 
D, = 0 Dar|r| = 0 (23) 
If we consider: 
a=0152m L=1,2m p = 1020 kg/m? 


So, after calculation, the hydrodynamic coefficients specific to this ROV will be grouped in the Table 1 and 
we will become: 


Table 1. ROVs hydrodynamic coefficients [S.I] 
D, 1401.56 Dwlu] 37.01 
D, 23258 Dlo] 186.04 


Dy 0 Dny|w| 186.04 
D, 0 Dnp|p| 10.04 
Dy 0 Dnglal 10.04 
D, 0 Dlr] 0 


2.3.4. External disturbances 

External disturbances affecting a marine vehicle are the wind (only for surface vehicles), waves (for 
surface or subsurface vehicles), and current, but the most important is the umbilical [17]. But in the presence 
of the embilical, the dynamic study should be more detailed [19]. Since it has no cable for our ROV and if we 
consider that the submarine will evolve in calm waters without the presence of waves or currents, these 
disturbances can be neglected. 
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2.3.5. Actuator effects 

They designate generally any power source which can generate a thrust exerted on the vehicle, in 
our case, they are four electric motors of brushless type mounted in configuration "x". Modeling the thrust of 
a thruster is relatively delicate because it depends on several parameters and is further complicated by the 
coupling of several thruster. Because in our case, the thrusts of the four thrusters are parallel to the axis Ox. 
seen from the rear, taking into account the configuration of the thrusters (Figure 5). The vector of forces and 
torques applied to the vehicle by the actuators is generally defined as [17]: 


T= (fo fy fale Tys T2)” (24) 
In our case, since we have four thrusters directed along the Ox axis, then: 

T = (fæ 0,0,0, Ty, ry (25) 
Therefore, the matrix of actuators: 


T=T.f; (26) 


Who becomes: 


1 1 1 1 
0 0 0 0 fi 
7 0 0 0 0 f 
T= 7 0 i i le (27) 
-d—h/2 -d—-h/2 +d—h/2 +d—h/2| lf, 
+d —d —d +d 


2.4. Simulation 
2.4.1. Numerical resolution 

The MATLAB language will be used for the simulation of the movement of this ROV, for that we 
will implement the program in this language. Similarly, and depending on the engine controls and also the 
response of the vehicle, the described trajectory will be visualized on a 3D interface. The parameters specific 
to the ROV of this study, as well as the equations which result from the kinematic, and dynamic studies will 
be integrated into the global program to visualize the various phases of the vehicle simulation. The ROV 
equation of state system can be defined as: 


nt) = Je). v(t) 


f = (28) 
te = M"1(T. f, — W (m) — D). v(t)) 
with: 
0 0 
Je(m) 0 0 0 
= 0 0 0 
0 0 0 Jalm) 
0 0 0 
and: 
X +m 0 0 0 0 0 
0 Y, +m 0 0 0 0 
we 0 0 Zy +m 0 0 0 a 
=-| o 0 0. Kotle 0 0 (30) 
0 0 0 0 M; + Ly 0 
0 0 0 0 0 N; +17 


Some numerical data specific to the ROV: 
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Jay = Jaz = 1p m.r? = 3,504 Kg. m? 
Jax = 0,143 Kg. m? 
rı = 60mm 


Tr, = 60mm (31) 
1=1,2m 

m = 28,6 Kg 

d=90mm 

h=15mm 


In the following part, we will use the numerical resolution of the differential equations to visualize 
the dynamic behavior specific to our submarine robot, of which the input parameters of the system will be 
only the different power configurations of the thrusters. The (28), which govern the model resulting from the 
kinematic and dynamic studies of the ROV, can be solved by several methods of numerical resolution of 
differential equations, but in this chapter, we have opted for the 4th order runge-kutta (RK) method for its 
remarkable precision in solving this type of equation [20]. 


2.4.2. Implementation 

This method is an efficient and reliable method for solving differential equations, and also can be 
used in a general way usable at any order of differential equations. To do this, it suffices to convert the order 
‘n’ of the equation to be solved to the number ‘n’ of equations to be solved. Nevertheless, it requires a 
meticulous sequencing for its application, because it involves four times more calculation than the Euler 
method for example. The sequence to apply this method for solving a Ist order differential equation with the 
4th order Runge-Kutta method, can be written as [21]: 


Y(t=0) = Yo 
Kı = f (to, Yo) -dt 
dt K 
K, = f (to +“, yo +) dt (32) 


dt K 

K3 = f(t +y) .dt 

K, = f (to + dt, Yo + K3) .dt 

yı = Yo +5 (Kı + 2K, + 2K; + K4) 


The diagram of this 4th order method above will be applied in our case for the solution of the 
equations. But we must first decompose our equation which is at the 2nd degree in 2 equations of the Ist 
degree. the coefficients of RK will be calculated simultaneously and alternately along the resolution process 
to finally find the desired values as: 


Y(t=0) = Yo h¢e=0) = ho (33) 
Kı = f (Yo) -dt Jı = g(ho) .dt (34) 
K = f (yo +4) dt J = 9 (ho +“) .at (35) 
K; = f (yo +Ë) dt Js = g (ho +“) .at (36) 
K, = f (Yo + J3) -dt Ja = g(ho + Kz) .dt (37) 
Yı = Yo +4. (Kı + 2K, + 2K; + Ky) hy = ho +4. U1 + 2)2 + 2J3 + J4) (38) 


Note: All these coefficients (vectors and matrices) are order 6. 

This algorithm will be implemented in MATLAB in the form of a program with the values of the 
thrusts as inputs. The most common combinations of thrusters will be tested to determine the reaction of the 
ROV in each case. 
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3. RESULTS AND DISCUSSION 

The system entries indicating the values of the thrusts developed by the thrusters will be given in the 
order [F;]-[F4]. To have the balance of the thrusts, normally they must have the same values. However, as the 
center of gravity of the vehicle is not on the axis (Ox), we will have an equilibrium coefficient of horizontal 
forces. So, this coefficient will indicate the proportionality between the horizontal thrusts and he is defined 
as: 


= D 
p=, Fa (39) 
We can see that if h = 0 we have a proportionality of p = 1 when F; = F,. For example: F} = 
22 and F, = 14 . To give a comparative idea between the study in an ideal and real environment, we will try 
to keep the same configuration and proportionality schemes previously entered. 


— Thrusters have proportionally the same thrust {18;18;28;28 } 

In this case (Figure 6) the robot keeps a straight path. Only the x position changes, the other five 
parameters remain zero. Its speed along the Ox axis increases initially, then stabilizes throughout the 
remainder of the simulation. 


Figure 6. Thrusters have proportionally the same thrust 


— Vertical thrusters have proportionally the same thrust {15;10;23;28} 
With this configuration (Figure 7), the ROV describes an arc in a plane parallel to the plane z = 0 at 
constant z, the radius depends on the power of the thrusters and only the angle theta varies. 


Figure 7. Vertical thrusters have proportionally the same thrust 
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— One of the four thrusters is proportionally different to the others {18;18;28;0} 
In this combination (Figure 8), the vehicle will describe a spiral whose radius and pitch will depends 
on the speed of the ROV and the deference of the power of this thruster with the other thrusters. 


Figure 8. One of the four thrusters is proportionally different to the others 


To confirm our results in practice, we started the realization phase in our Laboratory in the ENSEM 
Hassan II University. The production part is almost finished, prototyping and 3D printing are at a very 
advanced stage. In these figures (Figures 9 to 11), we have seen how the ROV is partly carried out, the 
assembly and the wiring are in an advanced stage, remain the dry tests, to then begin those in real conditions. 

Several methods are used to vary the signal of motor controls such as pulse width modulation 
(PWM) [21] (accessed May, 2021). 

(1] which is an efficient and widely used process including the naval sector, or with a control 
board dedicated to this task by using an URAT protocol (Universal Asynchronous 
Receiver/Transmitter) [22], [23]. , doi: 10.11591/ijpeds.v11.i4.pp1701-1710. 

However, we will first start with the control through the ESC (Electronic Speed Control) which is 
compatible with our brushless motors. For the control, we did not opt for an ambilical cable because it will 
limit the maneuvers of the robot [24], and we chose that the control of this robot will be remote and will be 
based on an arduino HC12 transceiver module. Its frequency range is between 433 MHz and 473 MHz which 
is already presented in a study [25], remains its adaptation for these modules for a greater depth. 


Figure 10. ROV’s assembly 


Figure 11. The control and wiring diagrams 
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4. CONCLUSION 

In this paper we were able to show theoretically that this submarine drone with four thrusters can 
move in a real aquatic environment without having to use rudders or diving bars. The various thrust 
combinations of the thrusters and the most common ones are tested and we have achieved the expected 
results. This justifies the choice of our solution. With the practical realization which is in progress, we can 
after confirm these theoretical obtained results and then proceed to the adjustments and finishes which are 
essential. In our next studies, we plan to expand the ROV control system with the tests and all the resulting 
optimizations. 
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